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Abstract. We investigate the strong asymptotics of Heine-Stieltjes polynomi- 
als - polynomial solutions of a second order differential equations with complex 
polynomial coefficients. The solution is given in terms of critical measures (sad- 
dle points of the weighted logarithmic energy on the plane), that are tightly 
related to quadratic differentials with closed trajectories on the plane. The 
paper is a continuation of the research initiated in [8]. tlowever, the starting 
point here is the WKB method, which allows to obtain the strong asymptotics. 



1. Generalized Lame equation 

In 1878 Heine studied the following problem: given two polynomials, 

p 

(1.1) A{z) = W{z — Qi) , B{z) = az^ + lower degree terms , a G C, 

i=0 

where the zeros of A are complex and pairwise distinct, describe the polynomial 
solutions of the generalized Lame differential equation (in algebraic form), 

(1.2) A{z) y"{z) + B{z) y' {z) - n(n + a - l)V,,{z) y{z) = 0, 

where Vn is a polynomial (in general, depending on n) of degree < p — 1; if deg V = 
p— 1, then V is monic. For p = 1 we can easily recognize in ( 1.2 1 the hypergeometric 
differential equation, while for p = 2 it is known as the Heun's equation (see |15j ). 
Heine [4] proved that for every n G N there exist at most 

+ p-l^ 
n 



(1.3) a{n) = 



different polynomials Vn such that ( 1.2 1 admits a polynomial solution y = Qn G P„; 
hereafter P„ stands for the set of all algebraic polynomials of degree < n. These 
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particular y„'s are called Van Vleck polynomials, and the corresponding polynomial 
solutions y = are known as Heine-Stieltjes (or simply Stieltjes) polynomials. 
Furthermore, if the polynomials A and B are algebraically independent then for 
any n G N there exist exactly a{n) Van Vleck polynomials V„, their degree is exactly 



p—1, and for each Vn equation ( 1.2 ) has a unique (up to a constant factor) solution 
y of degree n. According to recent results in jl6j . this statement still holds for 
arbitrary A and B for all sufficiently large n. 

Stieltjes discovered an electrostatic interpretation of zeros of the polynomials 
discussed in [4], which attracted common attention to the problem. He studied the 



problem ( 1.2 ) in a particular setting, assuming that all zeros of A are real and that 



all residues pk in 

B{x) _ Pk 



(14) 



A{x) a; - flfc 

are strictly positive. He proved in |18j (see also |20l Theorem 6.8]) that in this case 
for each rt G N there are exactly a{n) different Van Vleck polynomials of degree 
p—1 and the same number of corresponding Heine-Stieltjes polynomials y of degree 
n, given by all possible ways how the n zeros of y can be distributed in the p open 
intervals defined by the zeros of A. 

After these two contributions a vast number of research papers on this topic 
has been published, dealing mainly with the real situation considered by Stieltjes. 
Besides the relevant contribution to the algebraic theory |16j , we can mention here 
the paper f9], where the limit distribution of zeros of Heine-Stieltjes polynomials for 
the Stieltjes case has been established in terms of the traditional extremal problem 
for the weighted logarithmic energy on a compact set of the plane. 



A new approach to the asymptotics of the polynomial solutions of ( 1.2 1, based 
on a different type of equilibrium problem, has been developed in [8j , where zero 
distribution of Heine-Stieltjes polynomials was investigated. The present paper is 
in many senses a continuation of [8 . We use essentially the same tools: families of 
continuous critical measures and their representations in terms of quadratic differ- 
entials with closed trajectories. However, the starting point of our considerations 
here is different: instead of electrostatics we use the WKB method, which allows 
to obtain the strong asymptotics. 



2. Asymptotic formula for Heine-Stieltjes polynomials 

In this section we describe briefly the strong asymptotics of polynomial solu- 



tions of (1.2 1. As we will see immediately, it is based on the asymptotics of Van 
Vleck polynomials, so formally we should have started with the latter topic. We 
choose another way of presentation, which follows rather the logic of the proof of 
the main result. We formulate next a theorem on asymptotics of Heine-Stieltjes 
polynomials, which contains Van Vleck polynomials as unknown ("access") param- 
eters. 

Hereafter we assume that n G N is sufficiently large, so that there are exactly 
(T(n) Van Vleck polynomials, each of degree exactly p—1, see denote by Vn 
one of these Van Vleck polynomials, and let V„ stand for the set of zeros of Vn, 
A — {ao, . . . , ttp} the set of zeros of A, and ri„ C \ (^ U V„). 
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Fix any zq G C; in any simply-connected domain I? C f2 we can select a single- 
valued branch of ^JV^JA and define the natural parameter 

(2.1) ( = U^)= fxl^dt. 



In what follows, many statements will be made in terms of trajectories of the 
quadratic differential 

Vn{z) .2 

which are basically level lines of the function Im^„. 
We define also 




(2.2) Cn{z) = exp(^„(z)) = exp 
as well as the parameter 

(2.3) \n = n+ ^^y^ . 

Theorem 2.1. With the assumptions and notations above, let y — Qn be a Heine- 
Stieltjes polynomial (solution of ( |1.2| ) corresponding to Vn- Then: 

(i) there exist a set r„ = '^k>ilk.n comprised of a finite number of at most 

2p — 1 disjoint Jordan arcs jk.n, such that for an appropriate selection of 

the branch in ( |2.1| , 

(2.4) 2;(z) = ff„(z)C^(z) (l + £„(z)) , zeC\r„, 
where 

A , ( r B 



(2.5) F„(z) = ^— {z) j exp j (t) dt 

Moreover, en{z) — + with rt — > oo uniformly for dist(2;,r„) > C > 0. 
(ii) The endpoints of each jk^n belong to AU V„. Two arcs ^j_n and ^k.n for 
i j are either disjoint or have a common endpoint. The domain C \ r„ 
is connected, and y^Vn/A has a single-valued branch. Consequently, the 
branch of the square root in (2.4 1 is fully determined by condition 



lim zxl l'j^^^ = 1. 
A{z) 



Z — ^C30 



(iii) Each "fk,n can be selected as a part of a trajectory of the quadratic differ- 
ential 

(2.6) ^„ = -^(d,)^ 

which is close to its critical trajectory completed by two small segments 
connecting the trajectory with two points from AUVn- 



Observe that formula (2.4) is given in terms of polynomials Vn. In this way, 
the problem of asymptotics of Heine-Stieltjes polynomials is reduced to that of Van 
Vleck polynomials. We discuss the location of the zeros of VnS in Section |4j 
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The critical role is played by the set r„ = ^\^ilk,n- More exactly, formula (2.4) 
should be combined with the local asymptotic formula for in a neighborhood of 
7fc.n (but away from the endpoints): 

(2.7) Q„(z) = i/„(z) [C-{z) (1 + SnAz)) + C^"(2) (1 + e„.2(^))] , 

for certain small enj{z). Local formulas of this type can be obtained directly 
from the WKB estimates as explained in However, the existence of a global 

formula with the specified dominant term is not totally trivial. An arc 'yk,n enters 



the WKB analysis naturally as the Stokes line: the dominant term in (2.7) changes 
across jk.n (Stokes' phenomenon). This means that two terms are asymptotically 
equal on 7^. „ , and equation 

e"(^)--i 

defines the zeros of Q„. In other words, arcs jf^ n asymptotically carry zeros of (5„. 

The proof of the theorem is based on the LiouvuUe-Green (also known as the 
WKB) asymptotic formula. As we have mentioned above, the WKB method, the 
way it is explained in jllj . is ready to use in order to obtain local estimates. 
Construction of a global formula from the local ones constitutes a challenge, as 
it usually happens in the WKB analysis. So, the proof of the theorem, which is 
outlined next, is not entirely on the surface, and relies heavily on the fact that both 



the coefficients and a solution of (1.2) are polynomials 



3. WKB analysis of the Heine-Stieltjes polynomials 

The Liouville-Green (also known as the WKB) method is a standard tool of 
asymptotic analysis of ordinary differential equations; we refer the reader to the 
comprehensive account contained in the monographs [T and [llj. 

It was probably Nuttall [10 who applied the WKB analysis to the study of the 
strong asymptotics of the generalized Jacobi polinomials via a differential equation 



of the form (1.2 1 



3.1. Liouville's transformation. Recall that at this stage we assume that 
the Van Vleck polynomial Vn is fixed, and denote by y = Qn a corresponding 
Heine-Stieltjes polynomial of degree n. The following procedure is well known, so 
we summarize it here very briefly. 



Rewriting equation (1.2) in terms of 



(3.1) u{z)^y{z)exp(^l §^{t) dt 
we eliminate the first order derivative and get u"{z) = Fn{z)u(z), with 

(3.2) F„(z) = (A^ - p') \{z) + (^)\z) + (^)\z) , 



where A„ has been defined in (2.3) and 

P = -^T- 



1 /2 

Fix any zq S C; the Liouville-Green transformation w(^) — {£,') u, with ^ = £,n{z) 



defined in (2.1 1, yields 



(3.3) ^ = (A^.+.% (6)^^(6, 
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where 



(3.4) .g„(^) - G„(z(0) , G'„ = — ( — 1 + ( — ) - - ) - 



A I f bV f B 



and 



A detailed discussion of this transformation in the complex plain and of the resulting 
WKB approximation can be found, e.g. in 



3.2. Quadratic differentials on the Riemann sphere. At this point let 
us recall briefly some notions from the theory of rational quadratic differentials 
on the Riemann sphere; see [13j and [19 for further details. Given the quadratic 



differential Wn defined in (2.6 1, a smooth curve 7 along which 



= -Vn{z)/A{z) [dzf > ^ ImCTi(2) = const 



is a horizontal arc of ti7„; the distinguished parameter ^„ was defined in (2.1 1. More 
precisely, if 7 given by a parametrization z(t), t e (a,/?), then 

-^{<t)) (5)">0: («,/?). 

A maximal horizontal arc is called a horizontal trajectory (or simply a trajectory) of 
Wn- Analogously, trajectories of — n7„ are called orthogonal or vertical trajectories 
of Wn ; along these curves 

Vn{z)IA{z) {dzf > ^ ReCn(z) = const . 

We can define a conformal invariant metric associated with the quadratic dif- 
ferential w, given by the length element \dS^n\ = I M-^h the n7„-length of 
a curve 7 is 




(observe that this definition differs by a normalization constant from the definition 
5.3 in |19| ). Furthermore, if D is a simply connected domain not containing singular 
points of vjn, we can introduce the ron-distance by 

disi{zi, Z2;w,D) = inf{||7||^„ : 2:1,^2 €7, ^ d D} . 

Trajectories and orthogonal trajectories are in fact geodesies (in the tUn-metric) 
connecting any two of its points (see |13l Thm. 8.4]). 

A simply connected domain D not containing points from ^ U V„ is called a 
zun-rectangle if it is delimited by two horizontal and two vertical arcs of vj„; in 
other words, if ^n(-D) is a (euclidean) rectangle [a,b] x [c, d], and D 1-^ ^n{D) is a 
one-to-one conformal mapping. We call the value d — c the Wn-height, b — a the 
•n7„-length, and 2{b — a + d—c) the Wn-perimeter of D. Obviously, these definitions 
are consistent with the freedom in the selection of the natural parameter ^„ . 
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3.3. Local asymptotics. Theorems on local asymptotics of solutions are 
well-known (see [1] and using the local estimates that appear in Ch. 

VI, Theorem 11.1] it follows that if D is a n7„-rectangle such that function 5„ de- 
fined in ( |3.4| is holomorphic in $,n{D) = [a, b] x [c, d], then the differential equation 
( 3.3 1 has in [a, b] x [c, d] two linearly independent holomorphic solutions Wj, j — 1,2, 
of the form 

(3.5) w,{0 = exp {(-l)^+iA„C} (1 + e,(0) , 
and such that 

(3.6) \e,m<exp\— f^\g,,{t)\\dt\]-l, ^ e [a,b] x [c,d\ . 




The integrals here are taken following the progressive paths, i.e. contours along 
which Re(^) is non-decreasing (for j = 1) or non-increasing (for j — 2). In the case 
of the rectangle we may take Si = a + ic and S2 = b+id, so that the whole rectangle 
is reachable by progressive paths. Taking advantage of the fact that the coefficients 
of the original equation ( |1.2[ ) are polynomials, we can estimate the total variation 

A/(i^) =max / |5„(t)| \dt\ , 

where 7 is any horizontal or vertical segment in [a,b] x [c,d]. This yields the 
following result: 

Proposition 3.1. With the assumptions above, let D be an Euclidean rectangle 
such that S,n{D) C D and ^„ can be continued holomorphically to D in such a way 
that ^„(ylU V„) n Z? = 0. Let dn = dist(9-D, ^„(-D)) be the Euclidean distance from 
S,n{D) to the boundary of D. Then we can replace the estimates (3.6 1 by 

(3.7) 

In a -njn-rectangle D we can select a single valued branch of the function iJ„ 
introduced in (2.5). Then a direct consequence of the proposition above is 



Corollary 3.2. Let D be a w-rectangle. Then a general solution of (1.2) in D has 
the form 

(3.8) y{z) = H^{z) [/^lC^(^) {I + ei{z)) + K2C^-{z) (1-I-£2(^))] , 
with C,n defined in (2.2 1. We have 

(3.9) <exp(^^-^^ -1, zeD, 



where d„ is the Euclidean distance defined in Proposition 3.1 



3.4. Global asymptotic formula away from zeros. The result above shows 
that if we stay away from the singularities AyjV„, we can control the errors in the 
WKB approximation uniformly. This motivates the following definition. For £ > 
and i e C and subset if C C we denote 

DS) = {ze€:\z-t\<e}, D,{K) \J D{t, e) , 

and 
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Let Z? be a n7„-rectangle in C \ Dn,e containing infinity. We can select the 
branch of Cn and the initial point zq in (2.2 1 in such a way that |Cn(^)| > 1 + (5 for 
a (5 > 0, and there is a constant k such that 

(3.10) Q„(z)-iJ„(z) [C^(z) (l + ei(z)) + «C'"W (1 + £2W)] • 

It follows that if we continue all functions analytically in the maximal domain 
containing D such that |Cn(-2)| > 1 + f^, then formula (3.101 still holds. 

Next, standard arguments show that in a domain D, if Qn 7^ for all sufficiently 



large n, then there exists a dominant term in (3. 



Lemma 3.3. Let D = {z ^ C : a < Re(^„) < b, c < Im(^„) < d} be a simply 
connected zun-rectangle in C\ Dn.e, with 

TT 



d - c> 



|A„| 



If fl does not contain zeros of a solution y of (1.2) then for any choice of the lower 
limit of integration and with an appropriate choice of the branch of the square root 
in \2.^ we have 



(3.11) y{z) = H^{z)C-{z) (l + £„(z)) . 

Moreover, there exists a constant AI = M{e,D), independent of n, such that for 
any ■ujn-rectangle D' C D, 

M 



(3.12) 



\en{z)\ < 



An 



ten'. 



Now let us observe what happens if we have a union of two 7i7„-rectangles, D^^^ 
and D'-^\ such that D^^) n D^^) has an interior point, zq. Assume first that we 
take this zq as the lower limit of integration in (2.2 1 and (2.5 1. From Lemma 3.3 it 
follows that if y ^ 0, 



Z £ 



J 



1,2, 



y(z) = «:(^)il„(z)e(^) (l + £«(^) 
for certain non zero constants k*-^-* and k*^^-*. Evaluating at z — zq we get 
(1 + eilHzo)) = K^'^ (1 + e^^Hz„)) , z e D^^^ n Z?^^) ^ 



so that, 



(2) (1) l+£n''(^o) d) /, , / 

- (2)7— = Ml + £niz)) 



\£n{z)\ < 



M' 



In the intersection (neighborhood of zq) both expressions for y should match, hence 
we have that in D^^) n D^'^\ 

e"(z) (l + eW(z)) =C„^^"(^) {l + en{z)) (l + e^H^)) • 
Since D^^) n D^^) 

is an open set, necessarily the same branch of C„ has been taken 
in both sides of the previous identity. This argument shows that Lemma |3.3| is valid 
in the union D*^^' U D^^^ eventually with a different constant in the right hand side 
of ( [3I2I ). 

Now, modifying the lower limit of integration only changes the normalization 
constant k in (3.111. That means that in the representation (3.111 we can choose 
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any lower bound zq G -D'^' U D''^\ as long as we understand the function as the 
analytic continuation along a path in Z?'^) U D'-^K 

The discussion above motivates the following definition: 

Definition 3.4. Let D^^\ j = 1, . . . , fc, be a finite set of lUn-rectangles with Wn- 
hcight greater than |A„| and bounded t:t7„-perimeters. If each intersection D^^^ fl 
D^^'^^\ j ^ 1, . . . , k — 1, has an interior point, then their union Uj^iD^^^^ is called 
a finite Wn- chain (see Figure 



Figure 1. A tn„-chain. 



Thus, we have proved the following result: 
Lemma 3.5. Let D d C \ Dn,e be a domain, 
and simply connected Wn-chain in D. If a solution y of (1.2) does not vanish in D 



id D = U^^^iD^-J) C Z? a finite 



then for any choice of the lower limit of integration in D, and with an appropriate 



choice of the branch of the square root in (2.2) there exists a constant C such that 



(3.11)-(3.12) holds in D. The branch of (^n is obtained by analytic continuation in 



D. 

In other words, asymptotic representation can be continued along the chains 
of ti7„-rectangles, as long as we stay away from the zeros of the solution and of the 
sets A and V„. 

3.5. Zeros of a solution of the diffterential equation. Assume now that 
2/ is a nontrivial solution of (1.2 1 and zq G C \ is a zero of y. Let be a 
maximal simply connected Ci7„-rectangle in C \ I3„,e containing and let 7 be the 
vertical trajectory in 7 passing through zg- By assumption, zq is a regular point of 
the quadratic differential n7„, so that 7 is well defined. 



In f2 the expression (3.8 1 is valid. In particular, 

K2 I' 



£l(z) 



(3.13) y(z) = ^ 2A„en(^) = log, 
Hence, from the assumption y{zi^) = it follows that 

K2 1 + £1(2:0) 



mod (27ri) . 



2A„C«(2o) = log 



Ki 1 + £2(20) 



mod (27ri) 
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Assume that zi is another point on 7 that satisfies the following condition: 



(3.14) 



An 



\dUt)\ = K\hizo,z,)\L eN. 



where 7(2:07 zi) is the arc of 7 joining zq and zi. 

Definition 3.6. Let zq G C \ and let 7 be the largest connected vertical arc 

in C \ Dn^e of the quadratic difl:erential zun passing thorough zq. The set 

UJ = [jujj 



is a necklace in C\Dn^e corresponding to zq if all LUj, called beads, are tUn-rectangles 
of the form 

LUj^{zen:\ ReiUz) - U^M < <5, | IMU^) - U^M < S} , 
where each Zj satisfies condition (3.14), there exists a constant M such that X^6 < 
M, and all coj C C \ Dn,e- The vertical arc 7 is the string of the necklace. 



Figure 2. Necklace corresponding to zq. 
A direct consequence of Rouche's theorem is the following statement: 



Lemma 3.7. Let za £ C \ Dn^e be a zero of a nontrivial solution y of (1.2), and 
let Lu be a necklace in C\ corresponding to Zq. Then each bead Wj C C \ Dn^e 

contains one and only one zero of y. 

Furthermore, let be a simply connected Wn-rectangle in C\ Dn.e, and zq £ 



be a zero of a nontrivial solution y of (1.2). If to is the necklace in C \ D> 



corresponding to zq, then y{z) ^ for z £VL\uj. 

These results are applicable to any solution of ( 1.2 1. Now we concentrate on the 
n-th degree Heine-Stieltjes polynomial. The main difference is that we know that it 
has exactly n zeros, with account of multiplicity. This and Lemma |3 .7| immediately 
yield the following 

Proposition 3.8. Assume that zq G C \ Dn,t is a zero of a Heine-Stieltjes polyno- 
mial Qn, and let 7 be the largest connected vertical arc in C \ Dn^e of the quadratic 
differential vun passing thorough zq. Then 

II II ^ " + 2 

An 
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In particular, every string of a zero-carrying necklace has a finite Wn-length, and 
starts and ends at _D„ or is a closed curve. 



Indeed, by Lemma 3.7 every necklace carries at least A„||7||;^„ — 2 zeros of Qm 
where 7 is the string of the necklace. 

Summarizing, if Q„ is a Heine-Stieltjes polynomial, there exist a finite number 
of zero-carrying necklaces with the corresponding strings 7„j- such that all the zeros 
belong to the union of these necklaces or lie in 

4. Asymptotics of Van Vleck polynomials 



As it is clear from Theorem |2.1| the zeros of the Van Vleck polynomial Vn are 
the key parameters in the asymptotic expression for the Heine-Stieltjes polynomi- 
als. Hence, as a next step we derive a set of equation that will characterize their 
positions. We start with a formal argument and postpone to next section a more 
detailed discussion about consistency and meaning of these equations. 

Theorem 4.1. With the assumptions and notations above, let y — Q„ he a Heine- 



Stieltjes polynomial (solution of (1.2)) corresponding to Vn, and let ^n,k be the set 
of arcs defined in Theorem 2. 1 If there exists an e > such that all arcs 7fc,7i 
are disjoint and the Wn-distance between them is > s, then the following system of 
equations is satisfied: 



(4.1) — / . ——dt^ — -\ , fc=l,...,p-l, ruk&Mn, 

where the index se< M„ is a finite subset o/NU{0}. If for ^^^n we denote by rjk the 
set of endpoints ofjk,n that belong to A, then 

def sr^ B{a) A/ \ 

Furthermore, there exists a constant C — C{e) such that 

C 

\h,n\ < — ■ 

n 



This theorem is a simple consequence of the asymptotic formula (2.4), which 
is valid in a neighborhood of jk,n- Observing its increment along a closed Jordan 
curve 7fc_„ encircling 7^, „ in the positive direction and using the argument principle 
we get the formulas above. 



The system of equation (4.11 consists of p — 1 equations (since we have fixed 



the residue at infinity, equation k = p is dependent from the other p — I ones). 



More exactly, (4.1 1 presents a collection of systems of equations. In order to define 
it completely we need to specify: 

(a) the combinatorics: points from AUVn are arranged in p pairs (endpoints 

of 7fe,n)- 

(b) the homotopic types of curves 7fc_„ (once the combinatorics is fixed); and 

(c) the range of values the integer parameters in the right hand side of 



(4.1 ) may take. 



Observe that by (4.1 1 



1 f Vn{t) rnfe p,/2 / 1 , , ^ ^ 
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This system allows eventually to to find the position of "almost all" zeros of Van 
Vleck polynomials, with an error smaller than the distance between zeros. In most 
part of the range this error is 0(n"^). A complete analysis however is cumbersome 
and contains a combination of analytic, geometric and combinatorial arguments. 
We restrict our presentation here to the case p = 2 (three points), which at least 
has a trivial combinatorics and a rather simple geometry, with additional remarks 
on the case p — 3. 

4.1. Case oi p — 2. We have A = {qq, ai, 02} (in general, non-coUinear) and 
want to discuss the issues (a)-(c) raised above in order to define completely the set 
of equations determining the position of the zero of the Van Vleck polynomial. The 
first ingredient we need is the solution of the classical minimal capacity problem 
posed in the class of all continua in C containing A (the Chebotarev's problem). 
It was proved by Grotzsch [3] and Lavrentiev [6l \7\ that there exists a unique 
r* = r*(^) satisfying 

(4.2) cap(r*) = min{cap(-F) : F a continuum containing A}^ 

where cap(-) denotes the logarithmic capacity. If aj's are not coUinear, then T* is 
a union of three arcs, 70, 7i, and 72, connecting a point v* = v*{A) with points 
Oj's, respectively (see also [5j and [12j). We call this T* the Chebotarev's compact 
or Chebotarev's continuum corresponding to A, and point v* is the Chebotarev's 
center of the set A. If we define 

(4.3) M, = - / 
then 

Mo + Ml + M2 = 1. 
Observe that each Mj is the n7*-length of 7*, where 

^ = — (dz) . 

Next we define three analytic functions (elements) Wk{v), k = 0,1,2. We 
describe first wq as a, germ of an analytic function at v = qq, which allows unlimited 
analytic continuation to C \ For all w in a sufficiently small neighborhood of gq 

let Aq be the segment [ag, v], and Ai = [ai, a2]. Denote C \ (Ag U Ai). With 
R{z) == {z — v*)/A{z) we consider in fl the single- valued branch of ^/R given by 
the asymptotic condition limz^aa Zy^R{z) = 1, and define 

wo{v)^— [ ^/R(t)dt^^<f y^R(t)dt, 

where OAq is the doubly-connected component of dCl (the Caratheodory boundary 
of fi), with the boundary values of the branch of y/R specified above. It is clear 
that Wq is analytic at u = ao, and may be continued as a multi- valued analytic 
function to C Similarly, we define wi and W2, starting from v — ai and v = 02, 
respectively. The proof of the following statement can be found, e.g. in [5]: 

Proposition 4.2. For fc = 0, 1, 2 there exist an analytic arc ik, connecting Ok with 
V* , defined by 

4 = e C : Wkiv) e [0,Mfc]}. 



t - V* 



\dt\, /c = 0,1,2, 
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Furthermore, for k — 0, 1,2, Wfc(w) is univalent in a neighborhood Ck of £k (which 
is mapped by Wk onto a neighborhood of the interval [0, A/^,] j. 

Now we are ready to formulate a theorem about the asymptotic location of the 
zeros of Van Vleck polynomials for p — 2. 



Theorem 4.3. Let y = be a Heine-Stieltjes polynom,ial (solution of (1.2\) 
with A = {ao,ai,a2}, corresponding to Vn{z) = z — Vn- Assume that there exists 
an s > such that dist(w„,f*) > e (where v* is the Chebotarev's center). Then 
for all sufficiently large n, there exist an index k G {0, 1,2} and an integer value 
< vrik < [j^A/n] such that the following equation is satisfied: 

(4-4) Wk{Vn)^ ^ + ^ , ^ +^' 

and there exists a constant C — C{e) such that 

C 



\5kA < 



n 



Obviously, this statement does not cover, roughly speaking, en zeros out of 
n + I possible zeros of Van Vleck polynomials, but it provides an error estimate 
of order C{e)n~^. The exceptional set of "missed" Van Vleck zeros may be made 
smaller (up to a constant) for the price of relaxing the error estimate; also the 
technical details become more cumbersome. 

4.2. Case oi p — 3. Now we turn to the problems (a) and (b) related to 



the system (4.11, namely, we discuss the combinatorics and the homotopic type of 



curves 7fe,„ from Theorem 4.1 we leave the issue (c) of the range of values of ruk 
in the equations to the following Section. 

Here we illustrate the situation considering the points aj from A at vertices of 
a rectangle whose height (vertical size) is smaller then length (horizontal size), see 
Fig. [3] with the results of some numerical experiments, as well as Fig. [4] where the 
corresponding Chebotarev's continuum is depicted. 

In this case Vn{z) = {z — wi,„)(z — V2^n)', its zeros are determined completely 



by a system of two independent equations of the form (4.1 1. 

We claim that in this situation there are exactly 8 + 1 homotopically different 
groups of systems, see Figure [6] In eight of them we integrate along two arcs 7fc,„, 
each connecting a point from A with a zero of 1^. In the remaining case (Figure |6] 
bottom right) one of the curves will connect both zeros of Vn- 



For the case depicted in Figure |6] upper left, system (4.1 1 is written as 

Wl{ci,C2) 



W2(C1,C2) 



1 




Vn{t) 


dt = 


mk.i 


1 


■ni 


L \ 


J A{t) 


An 


An 


1 




IVnit) 


dt = 


mj.2 ^ 


1 


7ri 




/ A{t) 


An 


An 



B{ai) , Sk 



fc,n, 1 
An 



^(Q2) '^j,K,2 
2A'(a2) An ' 



where (771^^1,771^^2) is a pair of integers. We postpone the discussion about the 
possible range of variation of these constants until next section. In order to claim 
that 

C 

|4,nj| < — 
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Figure 3. Zeros of Heine-Stieltjes (small dots) and of the corre- 
sponding Van Vleck polynomials (fat dots) for Uj's at vertices of a 
rectangle (p = 3) . 



we need to leave aside agai n a fraction of zeros of Van Vleck polynomials (for which 
conditions of Theorem 2.1 are not satisfied), which is typically o(n^) (recall that the 



cardinality of the set of all Van Vleck polynomials corresponding to a Heine-Stieltjes 
polynomial of degree n is a{n) = 0{n'^)). 

5. The Van Vleck set (outline of the proof) 

This section is essentially based on the results of our recent paper [8], and we 
will be regularly referring the reader to some parts of this work for further details. 
We have seen that even in the simplest cases the geometry behind the set of 



equations (4.1 1 is quite involved; furthermore, we have not clarified yet the selection 
of the index set M„. However, the situation becomes much more clear when we 
take limit as n oo. 

We start by considering an arbitrary p for the price of omitting some non- 
essential details; a more thorough analysis will be carried out at the end for the 
particular (but non-trivial) case p — 2. 
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Figure 4. Chebotarev's compact corresponding to 4 points form- 
ing a rectangle. 



5.1. Zero distribution of Heine-Stieltjes polynomials and critical mea- 
sures. It is known (see e.g. [16 ) that the zeros of the Heine-Stieltjes polynomials 
accumulate on the convex hull of A. Hence, without loss of generality we may 
assume that there is a subsequence A C N such that the Van Vleck polynomials 
have a limit, 

p-i 

(5.1) liml4(z) = y(z) = rr(z-i;,)- 

This limit induces the following rational quadratic differential on the Riemann 
sphere, 

(.,2) "-W)^"''^ 

Let us consider the corresponding sequence of Heine-Stieltjes polynomials Qn, 
n G A; from Proposition |3.8| it follows that the zeros of Qn's lie asymptotically 
on critical trajectories of zu, and that these trajectories have the total n7-length 1. 
Further results can be obtained using the electrostatic interpretation of these zeros; 
it has been proved in [8 that any weak-* limit of the normalized zero-counting 
measure for 's is a continuous critical measure with respect to the set A of fixed 
points on the plane. Let us recall the definition in its simplest form, sufficient for 
what follows. 

With every (real-valued) Borel measure on C we associate its (continuous) 
logarithmic energy 

(5.3) E{fi) [ [ log ^^ d^i{x)d^l{y) . 

J J \x-y\ 

Any smooth complex-valued function h in the closure of a domain containing 
A generates a local variation of il by z i— s- z* = z -I- i h{z), t g C, and consequently, 
a variation of sets e i— > e* {z* : z e e}, and (signed) measures: /i i— > /i*, defined 
by /x*(e*) = ^(e). 
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Figure 5. Possible homotopic classes of curves jn.k for the case 
p — 3; compare with the examples depicted in Figure [3] 



We say that a signed measure fi supported in f2 is a continuous A- critical if for 
any h smooth in \ ^ such that = 0, 



(5.4) 



d 
dt 



t-^t) t 



^-critical measures appear quite frequently in many problems of approximation 
theory; any such a measure /i can be characterized in terms of its logarithmic 
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potential 

dcf 



as follows: 



Lemma 5.1 ([8 , Section 5.3). The logarithmic potential of an A-critical measure 
fi satisfies the following properties: 

(i) if supp(/i) ~ 7i U • • • U 7s, where 7j are the connected components of 
supp(/i), then 

U^{z) = Wj = const, 2^ € 7j, j = 1, . . . , s. 

(ii) at any regular point z G supp(/i) (that is, such that locally at z, supp(/i) 
is a smooth Jordan arc), 

(5-5) (^) = (^) ' 

where n± are the normal vectors to supp(/i) at z pointing in the opposite 
directions. 

Additionally, if z ^ supp(/i) \ ^ is not regular, then 
(5.6) gradC/^(z) = 0. 

Reciprocally, assume that a finite real measure fi, whose support supp(/i) consists of 
a union of a finite set of analytic arcs, supp(/i) = 71 U • • • U 7^, satisfies conditions 
(i) and (ii) above. Then /i is A-critical. 

In other words, (i) says that A-cv\i\ca.\ measures are in fact equilibrium mea- 
sures in a piece- wise constant external field, exhibiting additionally (see (ii)) the 
so-called S'-property, introduced first by Stahl |17] and in a more general context, 
by Gonchar and Rakhmanov [2] (where it was used to establish the well-known 
"1/9" -conjecture in approximation theory). A first rigorous proof of the connec- 
tion of a critical measure with the ^-property appeared in [14| . 

With a polynomial Q{z) = nfe=i(-^ ~ Cfe) of degree n we associate its (normal- 
ized) zero counting measure 

k=l 



The differential equation (1.2 1 is an expression of the fact that the zeros of Qn sit 



in equilibrium (zeros of the gradient of the discrete total energy) in presence of an 



external field depending from the residues pk in ( 1 .4 1 . The intensity of this field 



decays with n proportionally to 1 /n, so that in the limit we get 

Proposition 5.2 (0, Section 7). Let — ^{Qn) &e a zero-counting measure 
corresponding to a sequence of Heine- Stieltjes polynomials Qn. Then any weak-* 
limit point fi of is a unit continuous A-critical measure. 

In other words, weak-* limits of the normalized zero counting measures of QnS 
are unit positive ^-critical measures. The inverse inclusion (that any unit positive 
^-critical measure is a weak-* limit of the normalized zero counting measures of 
Heine-Stieltjes polynomials) is also valid, but it cannot be established within the 
framework of this paper. 
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Continuous critical measures can be also related with rational quadratic dif- 
ferentials on C by means of a variational argument. Our WKB analysis allows to 
establish this link, at least partially, more directly: 



Theorem 5.3. Let y = be a Heine- Stieltjes polynomial (solution of (1.2)) 



corresponding to Vn, and assume that (5.1) holds. Then the quadratic differential 



vj — — -(z) dz^ 

is closed (all its trajectories are either critical or closed), and there is a probability 
A-critical measure ii on C such that the normalized zero counting measures v{Qn) 
converge (in a weak-* sense and along the subset K) to ^. The support T — supp(/i) 
consists of critical trajectories of w, C \ F is connected, and we can fix the single 
valued branch of \JV j A there by lim^^oo ^\J^i,^^ l-^i^) = 1- With this convention, 

(5.7) lim|g„(z)|^/" =exp 




locally uniformly in <C\T , where a proper normalization of the integral in the right 
hand side is chosen, so that 

^(t)dt-log|z|j =0. 

The analysis above yields the following addendum, that we state using the 
notation and assumptions of Theorems |2.1| and |5.3[ 




Proposition 5.4. The support T ~ supp(/i) is comprised of p analytic arcs 7^, 
such that 

lim r„ = r. 

5.2. Structure of the family of positive critical measures. Theorem 
|5.3| is essential for understanding the asymptotics of Heine-Stieltjes polynomials, 
although this result is in a certain sense implicit, since it depends on the limit V 
of the Van Vleck polynomials Vn , that constitute therefore the main parameters of 
the problem. We must complement this description with the study of the set of all 
possible limits V. 

To this end it is convenient to consider a correspondence between closed qua- 
dratic differentials and .4-critical measures, independently of their origin. In gen- 
eral, this is not a one-to-one correspondence, since many critical measures may 
correspond to the same closed quadratic differential. It was pointed out in [8] that 



the bijection between closed quadratic differentials of the form (5.2) and signed 
.A-critical measures is restored if we restrict ourselves to signed measure with a 
connected complement of the support. Indeed, any such a measure is supported 
on a finite union of analytic arcs, and these arcs are necessarily critical trajectories 



of a quadratic differential like in ( |5.2[ ). Reciprocally, from (5.2) we can construct 
explicitly (using the Sokhotsky-Plemelj formulas) a measure fj, that satisfies (5.4). 

Thus, any closed quadratic differential uniquely generates an ^-critical mea- 
sure fi with a connected complement to its support. In general, such a is a signed 
measure, and for our purposes one has to select those quadratic differentials asso- 
ciated with positive measures. One of the ways to solve this problem goes through 
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describing the global structure of the trajectories of closed rational quadratic dif- 
ferentials with fixed denominators on the Riemann sphere, and the corresponding 
parameters (numerators), and later extracting from this set the subset giving rise 
to positive unit yl-critical measures. 

This way is implemented in for p — 2 at the end of this paper. For an arbitrary 
p such an investigation would be more difficult to carry out, in particular due to 
the heavy combinatorics involved. Below we present a direct construction from [8], 
based on the u-local coordinates that we introduce in the space of the y^-critical 
measures, as well as on the topological structure of the Chebotarev's compact for 
A. We start by introducing the w-coordinates. 

Let us recall the notation. We have the fixed set A = {ao, ai, . . . , a^} of distinct 
points on C, A(z) = n^=o(-^~ "^i)- ^'^^ signed ^-critical measure /i there exists 
a rational quadratic differential w on the Riemann sphere C given by 

(5.8) .^(z) = -i?,(z)(dz)^ ^W = ^. V{z)'='J[iz-v,), 

such that supp(/i) = = F = 71 U • • • U 7p is a union of trajectories of w. 

We begin by introducing the local coordinates under the assumption that mea- 
sures are in general position; this notion of genericity (as opposed to some more 
special or coincidental cases that are possible) means in our context that F is com- 
prised of exactly p disjoint arcs. In consequence, zeros Vj's of V are simple, and 
vector V {vi, . . . , Wp_i} G C*'^^ can be used as a local coordinate. Define also 

V = { V : ijj is closed}. 

Next we use these local coordinates to introduce the period mapping for our 
quadratic differentials. The Caratheodory boundary of C \ F consists of p compo- 
nents 7fe =^ U 7^7, with a positive orientation with respect to C \ F. We can 
consider 7fc as cycles in C \ F enclosing the endpoints of ■ Part of TZ over C \ F 
splits into two disjoint sheets, so we may consider 7^ as cycles on TZ. 

Let us define 

(5.9) Wk{v) = Wk{v,T)'^^ (£ y/R{z)dz, fc = l,...,p, 



where V R\~ are the boundary values of the branch of a/R in C \ F defined by 

lim^^oo ZyjR{z) — \. Clearly, the boundary values {\/R)± on 7^ are opposite 
in sign. Therefore, with any choice of orientation of 7^ and a proper choice of 
^/R = (\/i?)+ on 7fe, we will have 

(5.10) Wk{v)^— j ^/mzjdz, k=l,...,p. 



By the Cauchy residue theorem we have that wi + ■ ■ ■ + Wp = 1 for any v G 
(j~<p-i^ Thus, we can restrict the mapping v w to p — 1 components of w 
(wi, . . . , Wp-i) G C^^^. In this way, we have defined the mapping 

(5.11) V{-,T) : CP-^ ^CP~^ such that V{v,T) = w{v,T). 

Each component function Wj{vi, . . . ,fp_i) is analytic in each coordinate Vk (even 
if Vk is at one of the endpoints of 7^). Once defined by the integral in (5.101, this 
analytic germ allows an analytic continuation along any curve in C \ y^. Arcs 7^ 
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are not an obstacle for the continuation since the integral in (5.101 depends only 
on the homotopic class of F in C \ U t;). The homotopy of F is a continuous 
modification of all components simultaneously in such a way that they remain 
disjoint in all intermediate positions. Under this assumption we can continuously 
modify the selected branch of \/R in C \ F along with the motion of F. 

Proposition 5.5. Mapping w = 7-'(t>) is locally invertible at any v(vi, . . . G 
(C \ A)P-^ with V, ^ Vj for i^j. 

Furthermore, we have 

Proposition 5.6. Let jjP be an A-critical measure such that F = supp(/i'') has con- 
nected components . . . ,7^, and C \ F zs connected. Let zuq = RQ{z){dz)'^ be the 
quadratic differential associated with , where Rq = Vq/A, and v'^ ~ {vl, . . . , fp_i) 
is the vector of zeros 0/ Vq. Assume that is in general position (that is, all Vk 's 
are pairwise distinct and disjoint with A). Then for an e > and any rUj S K, 
j e {!,..., p~ I}, satisfying 

-/(7j)l < £, j = l,...,p-l, 
there exists a unique solution v Cz V of the system 
(5.12) Wj{v,T)=mj, j ^l,...,p-l. 

The quadratic differential vj = —R{z){dz)'^ , R{z) = 11^=1 (-^ ^ Vk)/A{z), is closed, 
and the associated A-critical measure 11 satisfies ^{'jj) = nij, j ~ 1, . . . ,p—l, where 
supp(/i) = 7i U • • • U 7p and supp(/i) is homotopic to supp(/x*'). 

This Proposition introduces a topology in the set of y^-critical measures in 
general position. We call a cell any connected component of this topological space. 
A measure in general position preserves sign along any connected component of its 
support. This sign is subsequently preserved when we homotopically modify the 
measure within its cell. In particular, if /i is a positive yl-critical measure, then all 
measures in the same cell with /i are positive. 

In ^-coordinates a cell G — G{T) in V is a subspace of ~ M'^p^-^, which is 
a manifold of the real dimension p — 1 , defined by p — 1 real equations of the form 



(5.13) luiwj{v)^lui\— I ^/R{i)dt\^Q, j = l,...,p-l. 




where F = 71 U • • • U 7p-i U 7p is a union of Jordan contours 7^ on TZ (double arcs) 
depending on v, but mutually homotopically equivalent for values of v from the 



same cell. Practically, any G G has a neighborhood of v satisfying (5.131 with 
constant 7fc's. 

Next, let V+ be the set of positive y^-critical measures; it will be comprised of 
3P~^ cells. Below we introduce these cells by means of the Chebotarev's continuum 
F* — F*(^), associated with A — {ao, . . . ,ap} (that is, the continuum solving 



the extremal problem (4.2)). It constitutes the "universal boundary", i.e., the 
intersection of all boundary points of all cells. In general, isolating a cell from its 
boundary point is a formidable task. We take advantage here that we are interested 
only in those that yield positive ^-critical measures, that are characterized by the 
trivial homotopy. 

The Chebotarev's continuum consists of critical trajectories of n7 = ~R*dz^, 
R* — V* /A. Assume that the set A satisfies the condition that A and V* do not 
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have common zeros and V*{z) ~ rifc=i(-^ ~ ^fc) does not have multiple zeros. Thus, 
the critical set A* — AU {vl, . . . , consists of 2p different points, and T* is 

comprised of 2p — 1 arcs, that are critical trajectories of va. Each trajectory joins 
two different points from A* . 

Now we come to the procedure of selection of combinatorial (rather than ho- 
motopic) types of cells; once the combinatorial type is fixed, the homotopic one 
will be determined from the Chebotarev's continuum, as described next. Each zero 
V* — vl, k — I, . . . ,p ~ 1, is connected by component arcs of F* with three other 
points, say al,a2,a^ G A*. We select one of these three arcs (for definiteness, 
[u*,a|]) and join two other arcs to make a single arc [02,03], bypassing v* (we 
think that the arc [02; '^3] still follows the two arcs from F*, but without touching 
V* , instead passing infinitely close to it). This procedure, carried out at each zeros 
vl of V* , creates a compact set F, and consequently, a cell G(F) of corresponding 
measures ^ G V+. 

The selection of F, and hence, of the cell G{T), is made by choosing one of the 
three connections for each v^; there are 3^""'^ ways to make the choice. Any choice 
splits F* into p "disjoint" arcs F* = 71 U • • • U 7p; out of them we select p — I arcs 
(to make an homology basis for C \F*) and then consider the corresponding cycles 
7fc, as described above. 

Finally, we describe the cell G(T) in terms of the mapping V. Let w = w(v) = 
7'(f,F); the cell G(F) is completely defined by the system 

(5.14) Wj{v) = ^lJ eR+, j = l,...,p-l; 

more precisely, there exists a domain A/(F) = {(/ii, . . . ,/Xp_i) G K^~^} such that 
for any point (^1, . . . , fJ-p-i) & M(F) system ( |5.14 l has a unique solution v € 



Moreover, the corresponding measure fi — iiy satisfies — i^j, and supp(/i) = 

71 U • • • U 7p = F^ is homotopic to F. 

Summarizing, a rough description of the set V+ of unit positive yi-critical mea- 
sures may be made as follows. The set V+ is a union of 3^^^ of closed bounded cells 
G(r) (F = 7iU- • •U7p may be selected in 3^"^ ways). The interior G(F) of each cell 
consists of measures fi in general position with supp(/i) homotopic to F. Interiors of 
different cells are disjoint. Chebotarev's measure fi* (Robin measure of F*) is the 
only common point of all boundaries: fi* = p|p9G(F). A graphical description of 
all these cells for p — 3, obtained from the construction just described, is contained 
in Figure [6] 

5.3. Case p — 2. In order to clarify the construction above let us discuss in 
more detail the simplest (but far from trivial) case of p = 2. 

Let us introduce the following set of the plane. For the quadratic differential 

(5.15) Wy = ^^--^ dz'^ 

A{z) 

define 

V {w G C : ti7^, is closed} , 
as well as the Van Vleck set 

V+ == {w G C : u is an accumulation point of the zeros of Van Vleck polynomials} . 



A direct consequence of Theorem |5.3| is that 

V+ C V. 
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As it was mentioned for the general case, this inclusion is proper. 

As in Subsection |4.1[ we denote by v* the Chebotarev's center of A (the value 



of V in (5.151 such that Wy* corresponds to the Chebotarev's continuum for A), 



and let A* ^AU{v*}. 

Theorem 5.7 (f8|, Section 8.3). The set V is a union of a countable number of 
analytic arcs £k, fc G each connecting v* and oo. 

Two arcs from V are either identical or have v* as the only finite common point. 
The homotopic type of the critical trajectories ofwy inC\A remains invariant on 
each arc £k \ A* . 

There are three distinguished arcs i^, k £ {0, 1,2}, such that 

(i) £k connects v* with infinity and passes through a^; 

(ii) for every v £ £k the homotopic class of trajectories of the closed quadratic 
differential vjy is trivial. 

In the terminology introduced in the previous subsection, this theorem says that 
each cell (connected component) in the topological space of ^-critical measures is 
homeomorphic to an analytic arc with endpoints either a,t v* , A or infinity. This 
settles the problem of identification of all closed quadratic differentials (fiy. Using the 
arguments described above we single out the differentials corresponding to positive 
measures: 

Theore m 5. 8. Let i^, k G {0,1,2}, be the distinguished arcs in V described in 
Theorem 5.7 The set V4. is the union of the sub-arcs of each li^, k G {0,1,2}, 
connecting with the Chebotarev's center v* (and lying in the convex hull of A). 

Furthermore, if v € £k H V+, k € {0,1,2}, then there is a critical trajectory 
7(w) of fly connecting v with the pole ak and such that 

(5.16) < fiyijiv)) < Mk, 



where is defined in (4.3 1. In this case both critical trajectories of Wy that 
constitute the support of fiy are homotopic to a segment. The bijection /jyi^^v)) ^ v 
is a parametrization of the set Ik H V+ by points of the interval [0, rrik]. 

In other words, the only three cells corresponding to positive y^-critical mea- 
sures are homeomorphic to three analytic arcs, joining the Chebotarev's center v* 
with each pole aj e A (in fact, these arcs form a star homeotopic and quite close 
to the Chebotarev's continuum, although in general not coincident with it). If v 
travels such an arc, the boundary of the cell is reached when either endpoint of 
this arc is met. If we continue further along the same arc, a new cell is entered, 
corresponding to sign-changing ^-critical measures (see an illustration of the cor- 
respondence between the position of t; on V and the trajectories of vjy in Figure 



Furthermore, by Theorem 5.8 the limit zero v oi V{z) = z — v satisfies an 
equation of the form 

1 r Im 



which is consistent with our construction in Subsection |4.1[ We can go back now 



to the set of equations (4.1) and add the following result: 
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Theorem 5.9. For any e > 0, equations 

(5.17) wk{v,.k) = ^ + ^, j=0,...,Mfe = [mfcn(l-£)], 

with fc = 0, 1, 2, uniquely define n = (1 — e)n± 1 points Vj^k- There exists a constant 
C = C(e) > such that any Van Vleck zero lies in a Cjr? distance from a point 

In this way, each zero of Vn "belongs" to a 0{n^'^) neighborhood of one of the 
points defined by the equations above. We need to show additionally that (again, 
up to a small neighborhood of the Chebotarev's center v*) all neighborhoods of 
this form contain at least one Van Vleck zero (establishing in this form a bijection 
between points Vj_k and the set V„). However, the techniques developed in this 
paper do not allow to complete this proof, so we formulate it as an open question: 

Conjecture 5.10. There exists a constant C — C{e) > such that to each point 
Wjjj defined by equations (5.17) it corresponds at least one zero w of a Van Vleck 



polynomial such that \v — Vj^k\ < C/n . 
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Figure 6. Position of v on £oU^iU^2 (left) and the corresponding 
trajectories of the differential voy in (|5.15[). 



